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Abstract 

The next-to-leading-order (NLO) cross section for the production of heavy quarks at large 
transverse momenta in 77 collisions is calculated with perturbative fragmentation functions 
(PFF's). This approach allows for a resummation of terms cx a s ln(p^/m 2 ) which arise in NLO 
from collinear emission of gluons by heavy quarks at large p± or from almost collinear branching 
of photons or gluons into heavy-quark pairs. We present single-inclusive distributions in p± and 
rapidity including direct and resolved photons for 77 production of heavy quarks at e + e~ colliders 
and at high-energy 77 colliders. The results are compared with the fixed-order calculation for 
m finite including QCD radiative corrections. The two approaches differ in the definitions and 
relative contributions of the direct and resolved terms, but essentially agree in their sum. The 
resummation of the a s ln(|j^/m 2 ) terms in the PFF approach leads to a softer p± distribution 
and to a reduced sensitivity to the choice of the renormalization and factorization scales. 
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1. Introduction 



A large number of equivalent photons is generated at high-energy e + e colliders, giving rise to 
the production of heavy-quark pairs in two-photon collisions. This process has been studied for 
charmed particles by several experiments at PETRA, PEP, TRISTAN, and LEP [gj. At LEP2, a 
total of ~ 350 000 cc and ~ 1500 bb pairs will be produced in 77 collisions for an integrated luminosity 
of / C = 500 pb _1 [0]. The yield of heavy quarks at high-energy e + e~ linear colliders is even higher, 
depending in detail on the spectrum of the beamstrahlung photons, which strongly varies with the 
machine design and operation. If the novel method of Compton back-scattering of laser light can 
be made work Q, it will be possible to generate high-luminosity beams of real photons carrying 
~ 80% of the electron/positron energy. The high-statistics data from the future experimental 
facilities will allow for a detailed comparison of the next-to-leading order (NLO) predictions with 
experimental results not only for total production rates, but also for various differential distributions. 
These analyses will provide us with information on the dynamics of heavy-flavour production in a 
kinematical range very different from that available in 77 collisions at present colliders. 

Three mechanisms contribute to the production of heavy quarks in 77 collisions: (i) In the 
case of direct (DD) production, the two photons couple directly to the heavy quarks. No spectator 
particles travel along the 7 axes, [it) If one of the photons first splits into a flux of light quarks 
and gluons [|J], one of the gluons may fuse with the second photon to form the QQ pair. The 
remaining light quarks and gluons build up a spectator jet in the split-7 direction (single-resolved 
(DR) 7 contribution). The total 77 cross section of this mechanism depends on the parton density 
functions (PDF's) of the photon (Hi) If both photons split into light quarks and gluons, the 

QQ pair is accompanied by two spectator jets (double-resolved (RR) 7 contribution). Since the 
photon PDF's scale as aaj 1 , the DR and RR processes are of the same order as the DD process. 

Many features of the above-mentioned production mechanisms are calculable in perturbative 
QCD. The mass of the heavy quark, m 3> AqcDj ac ts as a cutoff and sets the scale for the per- 
turbative calculations. The production cross section factorizes into a partonic hard-scattering cross 
section multiplied by light-quark and gluon PDF's M. Inherent in this factorization scheme is the 
notion that the only quarks in the photon are the light ones. There are no contributing subpro- 
cesses initiated by an intrinsic heavy flavour coming directly from the photon PDF's. In leading 
order (LO), direct production is described by the partonic reaction 7 + 7 — > c + c while the resolved 
contributions involve the channels g + 7 — > c + c (DR) and q + q — > c + c, g + g — > c + c (RR) , where 
q are light (massless) flavours ||. The NLO corrections to these processes have been calculated and 
found to be substantial @,[H],|ll]@ • A comparison between the NLO results and experimental data 
on the total cross section of charm-quark production in two-photon collisions has shown satisfactory 
agreement 13]. 

One might expect that the massive approach is reasonable only in those kinematical regions 
where the mass m and any other characteristic energy scale like p± are approximately of the same 
magnitude and significantly larger than Aqcd- In NLO, terms oc a s ln(pj_/m 2 ) arise from collinear 
emission of gluons by heavy quarks at large transverse momenta (p±) or from almost collinear 
branching of gluons or photons into heavy-quark pairs. These terms are not expected to affect the 
total production rates, but they might spoil the convergence of the perturbation series and cause 
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large scale dependences of the NLO result at p± 3> In the massive approach, the prediction 
of differential cross sections is thus limited to a rather small range of p±_ ~ m. An alternative 
way of making predictions at large p±_ is to treat the heavy quarks as massless partons. The mass 
singularities of the form \n(p\/m?) are then absorbed into the PDF's and fragmentation functions 
(FF's) in the same way as for the light u, d, and s quarks. The crucial difference to the production of 
light hadrons is the fact that the initial-state conditions for the heavy-quark FF's can be calculated 
within perturbative QCD and do not have to be taken from experiment |jq[. Such perturbative 



fragmentation functions (PFF's) have been been used in Ref. [16] to study the production of large-pj_ 
bottom quarks in pp collisions. Meanwhile, similar analyses have been carried out for charm-quark 
production in photon-proton collisions at HERA [17,18|. In the massless scheme, the heavy quark is, 
of course, considered to be one of the massless active flavours in the photon PDF's. We expect the 
massless PFF approach to be better suited for the calculation of the differential p± distributions at 
NLO in the region p± 3> m. The smail-pj_ region is, however, not calculable without retaining the 
full m dependence. The massless cross section diverges in the limit p± — > 0, and total production 
rates can not be predicted. 

At this point, it is worth mentioning that the massive and the massless approaches are expected 
to give different descriptions of the DD, DR, and RR contributions. In the massive approach, 
the direct and resolved contributions are separately well defined through NLO. In fact, there is 
no need to perform any factorization for the incoming photons. This means that an experiment 
analyzing data without hadronic activity in the directions of the incoming photons would directly 
probe the NLO calculation of DD heavy-quark production. This is, however, not true in the massless 
approach. If the charm is treated as a massless parton, the photon splitting to cc must also undergo 
a subtraction procedure. The DD piece then becomes factorization-scheme dependent and looses 
its direct physical interpretation. Consequently, the same is true for the DR and RR parts. As 
a matter of principle, the three contributions cannot be experimentally separated any more with 
NLO accuracy. Only their sum corresponds to a physical observable. On the other hand, if one 
were to extend the massive approach up to large p±, the DD prediction would, of course, still be 
unambiguously defined, but it would be affected by large logarithmic terms, which would render 
it unreliable. Only at moderate p± ~ m and through NLO, it is therefore possible to probe the 
theoretical predictions for direct and resolved contributions separately by an experimental analysis. 

In the near future, we expect experimental data in the intermediate pj_ range, where p± > m 
rather than p± 3> m. Then, the problem how to proceed in this p± region arises. In order to 
investigate the region where p± > m, we calculate the differential cross section d 2 a/dydp\ as 
a function of p± with fixed rapidity y. We compare the results in the two approaches: (i) the 
massive-charm approach with m = 1.5 GeV, in which we compute the cross section for open charm 
production, and (ii) the massless approach, where we evaluate the same differential cross section 
for inclusive charm production using PFF's. In both calculations, we include the DD, DR, and RR 
processes up to NLO. The massive calculation is based on the work presented in Q, in which the 
NLO theory for DD and DR production was elaborated. The calculation of the massive RR cross 
section relies on the work of |10jiT|. The massless calculation proceeds along the lines of |16-|is|1 on 
the basis of the DR and RR hard-scattering cross sections obtained in [20,2l|]. The NLO corrections 

1 Similar potentially large terms oc a s ln(Q 2 /m 2 ), Q being the photon virtuality, do appear in NLO calculations of 
heavy quark electroproduction cross sections. The resummation of these logarithms is being considered in a series of 
papers 
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to the DD cross section were derived in |2^] and recently confirmed in |23|]. In the massless approach, 
we adopt the PFF's calculated in [15]. With this choice of PFF's, the LO results in the massive 
scheme approach the LO massless results in the limit m — > 0, if we restrict ourselves to the same 
parton subprocesses. In LO, the PFF's are equal to 5(1 — x) (see Eq. ([!])) showing that, in LO, we 
have full correspondence between the massless and massive approaches. However, they must differ 
in NLO, where the limit m — > is not possible due to the unabsorbed mass-singular ln(pj_/m 2 ) 
terms in the massive approach. In the massless approach, these terms are contained in the higher- 
order terms of the photon PDF's and the PFF's. This situation is very similar to our previous 
study for photon-proton collisions at HERA [17.18]. Indeed the direct and resolved contributions in 
■yp scattering correspond to the DR and RR contributions in 77 collisions with the proton PDF's 
replaced by the photon PDF's. In addition, we now have the DD contribution. 

The outline of our work is as follows. In Section 2, we shortly describe the basic formalism for the 
massless PFF approach. Section 3 contains the numerical results. Our conclusions are summarized 
in Section 4. 



2. The perturbative fragmentation function approach 



In this section, we describe in some detail the PFF approach to heavy-quark production at large 
p±. This technique was proposed in Jllj J and first applied in the context of hadron collisions for 
describing large-pj_ bottom production The basic assumption is that when a large scale is 

governing the production process (in our case, this is p±, with p±_ 3> m) the heavy quark itself is 
produced as if it was massless. Technically speaking, non-singular mass terms are suppressed in the 
cross section by powers of m/p±. The important mass terms appear whenever the virtuality of the 
heavy quark is small. This happens in the initial state when the heavy parton is emitted from the 
colliding hadron, and in the final state when the partons materialize into a massive quark. 

This qualitative picture of heavy-quark production in the large-pj_ limit can be substantiated 
at NLO in the following way: (i) The hard-scattering cross sections are calculated in the massless 
approximation, and the collinear singularities are subtracted according to some factorization scheme, 
e.g., the MS scheme. Since the heavy parton is taken to be massless, also the singularities arising 
from its splittings are subtracted, (ii) As for the initial state, the heavy parton is accommodated 
in the PDF's like a light flavour. The massiveness of the quark is usually taken into account by 
including it in the evolution only above a scale set by its mass. (Hi) As for the final state, the PFF's 
characterize the hadronization of the massless partons into the heavy-quark state. Exploiting the 
fact that the produced quark has mass m 3> Aqcdi universal starting conditions for these FF's can 
be calculated within perturbative QCD (therefore they are denoted PFF's) at a scale fiQ of order 



m. In [15], these starting conditions were calculated at NLO in the MS scheme. They read 

a s (fi ) ^ {l + x 2 



D Q (x,no) 



5(1 -x) + 
a s (/J,o 



2tt 



-C F 
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7^ + (I-*) 2 ] log 4, 



0, 



(1) 



where D® refers to the fragmentation of parton a into the heavy quark Q, Cp = 4/3, and Tf = 1/2. 
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(iv) Finally, the PDF's and PFF's are evolved in NLO up to the chosen factorization scale (which is 
usually of the order of p±) via the Altarelli-Parisi equations and convoluted with the hard-scattering 
cross sections. Notice that, in this approach, the heavy-quark mass enters the calculation only via 
the starting conditions of the PDF's and PFF's. 

In this framework, the large logarithmic terms are resummed in the following way. The would-be 
mass singularities oc ln(p^_/m 2 ) are split into two parts. One part oc \ri(p\/n 2 ), where [i is the 
factorization scale, appears in the hard-scattering cross sections, which have no dependence on m. 
This part may be eliminated by choosing \x ~ pj_. The other part oc In (\j? jm 2 ) is absorbed into the 
PFF's. The large In (fJ 2 /^), with fio ~ m and \x ~ p±, is implemented via the evolution equations, 
and therefore these large logarithms are resummed. The residual terms oc In (m 2 //xg) connected 
with the starting condition in Eq. ([l]) are treated at fixed order in perturbation theory. 



3. Results 

In the following, we collect our results for the cross section d 2 a/dydp 2 L for three cases of partic- 
ular interest: (i) LEP2 with y/s = 175 GeV, (ii) Next Linear Collider (NLC) with y/s = 500 GeV 
assuming the TESLA design, and (Hi) NLC with y/s = 500 GeV operated in the 77 mode imple- 
mented by backscattering of laser light on the e + and e~ beams. 

We start with our analysis relevant for LEP2. The quasi-real-photon spectrum is described in 



the Weizsacker- Williams approximation (WWA) by the formula [24 



_ a f l + (l-x) 2 E 2 9 2 (l-x) 2 + m 2 e x 2 
h[X) ~27r) x m m 2 x 2 +{ 



TT? e X 1 



E 2 9 2 (l — x) 2 + m 2 x 2 x 



(2) 



where x = E^/E e is varied over the full range allowed by kinematics and 9 C is the maximum 
angle under which the outgoing electrons (positrons) are tagged. In our LEP2 analysis, we choose 
9 C = 30 mrad. All calculations are performed at NLO in the MS renormalization and factorization 
scheme using the two- loop formula for a s . The DR and RR cross sections are calculated using the 
photon PDF's of Gliick, Reya, and Vogt || transformed to the MS scheme. The renormalization 



and factorization scales are set to \x = J p\ + m 2 . It is clear that in the massive scheme only three 
flavours are active in the initial state and in the evaluation of a s , whereas in the massless scheme also 
the charm distribution in the photon contributes to charm quark production, and a s is calculated 
using four active flavours. In the massive calculation and in the PFF's, we set m = 1.5 GeV. 

In Figs, la-d, we show the DD, DR, and RR contributions to the cross section d 2 a/dydp 2 L as 
a function of p± for rapidity y = and their sum, respectively, both in the massless and massive 
schemes. However, one should bear in mind that, beyond LO, the separation into the DD, DR, and 
RR channels depends on the factorization scheme and scale and has no direct physical meaning. 
[As we have discussed above the DD channel in the massive scheme carries unambigous physical 
meaning.] Nevertheless, we consider these channels separately in order to assess their relative 
importance. As may be seen from these figures, the cross section in the considered p± range is 
dominated by the DD contribution. For the DD and DR components, the massless and massive 
calculations are very similar, whereas in the RR case the massive prediction is at least one order of 
magnitude smaller than the massless one. Comparing the two schemes, we see that the difference 
is somewhat larger for the DD contribution than for the DR one. In both cases, the massive cross 
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section exceeds the massless one, which is still the case for the total sum. The discrepancy increases 
for increasing pj_. The situation for the DR and RR contributions is comparable to what we have 
observed for photoproduction of charm quarks at HERA []l7| , |l8| . However, due to the DD channel, 
which does not exist in photoproduction at HERA, the total 77 results in the massless and massive 
schemes exhibit a pattern somewhat different from [l7| , |l8| ]. In total, we observe that the massless 
prediction for the total sum is smaller than the massive one, in particular for large p±. We attribute 
this to the presence of higher-order leading-logarithmic terms in the PFF's incorporated in the 
massless approach, which soften the p± distribution of the heavy quark. 

We now repeat the analysis of Fig. 1 for the y spectrum at p± = 10 GeV and display the DD, 
DR, and RR contributions as well as the total sum in Fig. 2a-d, respectively. Comparing the results 
in the massive and the massless approaches, we observe a pattern similar to the p± distribution. For 
given p± , the kinematically allowed y range in the massless theory is larger than that in the massive 
theory. This does not show up in Figs. 2a, b, and d due the influence of the PFF's. Moreover, at 
the edges of phase space, where only soft-gluon emission is allowed, the PFF's are not reliable due 
to the missing resummation of Sudakov terms. In the region of large cross section, the shapes for 
the two approaches are very similar. 

Next, we consider the predictions for the NLC with TESLA architecture. It is well known 
p5| that, at the NLC, photons are produced not only by bremsstrahlung but also via synchrotron 
radiation emitted by one of the colliding bunches in the field of the opposing one. This phenomenon 
is called beamstrahlung. The details of the beamstrahlung spectrum crucially depend on the design 
and operation of the NLC. For certain NLC concepts, beamstrahlung can jeopardize the overall 
physics potential due to severe smearing and lowering of the available centre-of-mass energy. In our 
study, we select the TESLA design, where the unwanted effects of beamstrahlung are reduced to an 
almost unnoticeable level. We coherently superimpose the WWA and beamstrahlung spectra. We 
compute the WWA spectrum from Eq. (0) with 6 C = 175 mrad and the beamstrahlung spectrum 



from the expression given in [25], with parameters T e g = 0.039 and a z = 0.5 mm [26]. The p± 
distributions for y = and yfs = 500 GeV are shown in Figs. 3a-d, again for the DD, DR, and RR 
components, and their sum, respectively. Apart from an overall enhancement due to the increased 
WWA logarithm, the p± spectra exhibit features very similar to the LEP2 case, including the 
relation between the massless and massive approaches. The corresponding y spectra are shown in 
Figs. 4a-d. Due to the admixture of beamstrahlung, their shape differs from that of the pure WWA 
case in Figs. 2a-d. This difference is most pronounced in the DD contribution, shown in Fig. 4a, 
which also governs the shape of the total sum. 

To achieve the highest possible photon energies with large enough luminosity, it has been pro- 
posed to convert the NLC into a 77 collider via backscattering of the e + and e~ beams on high- 
energetic laser light ||. The corresponding photon spectrum is given by || 

, . , 1 / 1 4x Ax 2 \ 

/ 7 W = 7^ 1 - X + , + 



G(k) \ 1-x k(1-x) k 2 {1-x) 2 )' 

with 

Notice that this spectrum extends only up to x max = k/(1 + k). In order to avoid the production 
of e + e~ pairs in the collisions of the primary laser photons and the high-energetic back-scattered 
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photons, one needs to arrange the experimental set-up so that k <, 4.83. In our analysis, we choose 
k = 4.83, so that x max = 83%. In Figs. 5a-d, we show the p± distributions for y = due to the DD, 
DR, and RR channels, and the total sum, respectively. Compared to the previous cases, the relative 
magnitudes of the DD, DR, and RR components have changed. At small p±, around 5 GeV say, 
the massive cross section is dominated by the DR component, which makes up approximately 60% 
of the total sum, while the DD component is completely negligible. At the same p± values in the 
massless case, RR is the largest component, being 65% of the full result. The DD contribution is 
again negligible. The massless cross section is only slightly larger than the massive one. At large p±, 
at around 20 GeV, the situation is different. In the massless case, the DD, DR, and RR components 
all have the same order of magnitude. In the massive case, the RR contribution is small, and the 
total cross section is built up by the DR and DD contributions approximately in the ratio 3:1. 
Looking at Fig. 5d, we see that p± distributions in the two schemes almost coincide. However, the 
massless result falls off somewhat more strongly with pj_ increasing. As may be expected, in both 
schemes, the DD contributions are insignificant in the low-p_|_ range but become important for high 
p±. The corresponding y spectra for pj_ = 10 GeV are plotted in Figs. 6a-d. As far as the relative 
importance of the individual contributions in the two schemes is concerned, we recognize the same 
pattern as in Figs. 5a, b, and c, which refer to y = 0. The line shapes are completely different from 
those encountered in the TESLA case. The y spectrum of the total sum peaks near the phase-space 
boundaries. This is caused by the different photon spectrum, which is now peaked at the upper 
edge. By contrast, the numerically small DD contribution is almost y independent. Although the 
total sums in the massless and massive calculations have different decompositions, they nevertheless 
agree very well with each other over the full y range. 

Finally, we investigate the scale dependence of our results. To this end, we introduce a dimen- 
sionless scale parameter £ and set all scales equal to £ J m 2 + p\ . In Fig. 7, we plot, versus £, 
the massless and massive cross sections at y = and p± = 20 GeV for LEP2, TESLA, and Laser 
spectrum. For 0.5 < £ < 2, we observe minor scale variations, below 15%. As expected, the scale 
dependence is somewhat reduced in the massless cross section. We checked that the situation is 
similar for p± = 10 GeV. 



4. Conclusions 

In this paper, we compared two approaches, the massless and massive schemes, for calculating 
inclusive charm quark production in three different arrangements of 77 reactions. The cross sec- 
tions were computed at NLO, with PFF's included in the massless case. Comparing the massless 
and massive results in the DD, DR, and RR channels, we found essential differences, which are 
compensated to a large extent in their sums. 

In all three arrangements, the cross section at large p± was found to be somewhat larger in 
the massive case. This is in line with the findings of [19] for bottom hadroproduction and of 



[18] for photoproduction at HERA. In the latter case, too, it was observed that, in the massless 
approximation, the direct and resolved contributions, which in our case correspond to the DR and 
RR contributions, are comparable. On the other hand, the resolved contribution of the massive 
calculation was found to be suppressed ]jT7| , pLql , which also agrees with our 77 results. 

Concerning the small p± region, we notice that the total massless result usually overshoots the 
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massive one. This is due to the missing mass terms in the massless approach, which render it less 
reliable in this region. This effect is enhanced with respect to what was observed in [plj|j , for the 
following two reasons: In this work, we use massless kinematics, whereas, in JlS] , massive kinematics 
was employed. Moreover, in [18], the photon PDF set ACFGP-ho (mc) |7| was used. In contrast to 



GRV used here, in this set, the massive Bethe-Heitler formula is built in for the starting condition 
of the charm PDF of the photon. This had the effect that, especially in the low-p± region, the cross 
section was reduced. 

In conclusion, we have found that, in the region of common validity, the massive and the massless 
approaches yield comparable results and display similar scale dependences. In the large p± region, 
the massless approach predicts a lower cross section, due to the resummation of a s ln(p^_/m 2 ) terms, 
and leads to a reduced sensitivity to the choice of the renormalization and factorization scales. 
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FIGURE CAPTIONS 

Fig. la d Inclusive cross section d 2 a/dy &p\ for e + e~ — > e + e~c/c + X as a function of p± for 
yfs = 175 GeV (LEP2) and y = in the massless (solid lines) and massive (dashed lines) schemes: 
(a) DD, (b) DR, (c) RR, and (d) total sum. 

Fig. 2a-d Inclusive cross section d 2 a/dydp 2 ± for e + e~ — ► e + e~c/c + X as a function of y for 
yfs = 175 GeV (LEP2) and p± = 10 GeV in the massless (solid lines) and massive (dashed lines) 
schemes: (a) DD, (b) DR, (c) RR, and (d) total sum. 

Fig. 3a-d Same as Fig. la-d for sfs = 500 GeV (NLC with WWA plus beamstrahlung) . 
Fig. 4a-d Same as Fig. 2a-d for = 500 GeV (NLC with WWA plus beamstrahlung). 
Fig. 5a-d Same as Fig. la-d for ^fs = 500 GeV (NLC with laser spectrum). 
Fig. 6a d Same as Fig. 2a-d for y/s = 500 GeV (NLC with laser spectrum). 

Fig. 7a— c Scale dependence of d 2 a/dydp\ for e + e~ — ► e + e~c/c + X at y = and p± = 20 GeV 
in the massless (solid lines) and massive (dashed lines) schemes: (a) yfs = 175 GeV (LEP2), (b) 
yfe = 500 GeV (NLC with WWA plus beamstrahlung), and (c) ^ = 500 GeV (NLC with laser 
spectrum). The renormalization and factorization scales are identified and set equal to £,Jp]_ + rn 2 . 
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Fig. 2 
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Fig. 3 
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Fig. 4 
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Fig. 5 
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Fig. 6 
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